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LONG-TERM  GOALS 

Development  of  a  physical  model  of  high-frequency  acoustic  interaction  with  the  ocean  floor, 
including  penetration  through  and  reflection  from  smooth  and  rough  water/sediment  interfaces, 
scattering  from  the  interface  roughness  and  volume  heterogeneities  and  propagation  within  the 
sediment.  The  model  will  aid  in  the  detection  and  classification  of  buried  mines  and  improve  SONAR 
performance  in  shallow  water. 

OBJECTIVES 

1)  Determination  of  the  correct  physical  model  of  acoustic  propagation  through  ocean  sediments  and 
scattering  from  sediment  interfaces  through  the  analysis  of  in  situ  measurements. 

2)  Development  of  predictive  models  that  can  account  for  the  all  of  the  physical  processes  and 
variability  of  acoustic  propagation  and  scattering  in  ocean  environments  with  special  emphasis  on 
propagation  in  shallow  water  waveguides  and  scattering  from  ocean  sediments. 

3)  Development  of  the  new  experimental  techniques  to  measure  geo-acoustic  parameters  in  the  ocean. 

APPROACH 

1)  Finite  Element  Modeling  for  Range  Dependent  Waveguides :  Finite  element  modeling  is  applied  to 
a  waveguide  with  measured  range  dependent  geo-acoustic  parameters  in  order  to  calculate 
transmission  loss  and  reverberation.  Two  different  approaches  to  the  three-dimensional  problem 
are  taken:  axially  symmetric  models  and  longitudinally  invariant  modeling. 

2)  Analysis  of  Normal  Incidence  Bottom  Loss  Measurements  for  Range  Dependent  Geoacoustic 
Parameters:  Bottom  loss  data  from  5-30  kHz  were  collected  as  part  of  the  Target  and 
Reverberation  Experiment  2013  (TREX13).  These  data  were  analyzed  and  range  dependent 
geoacoustic  parameters  were  derived  for  the  TREX  reverberation  site  including  bottom  loss  and 
scattering.  The  data  were  compared  with  multibeam  sonar  data  taken  at  the  same  site. 

3)  Measurements  of  Bottom  Loss,  Sediment  Structure  and  Interface  Roughness  at  the  Glider  Sensors 
and  payloads  for  Tactical  characterization  of  the  Environment  (GLISTEN  15)  Experiment:  Normal 
incidence  bottom  loss  measurements  from  5-20  kHz  and  interface  roughness  data  were  collected 
using  an  ROV  platform  as  part  of  the  GLISTEN  experiment. 

WORK  COMPLETED 

Finite  Element  Modeling  for  Range  Dependent  Waveguides 
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A  finite  element  model  of  propagation  and  scattering  was  calculated  based  on  the  domain  shown  in 
Fig.  1  and  the  geoacoustic  parameters  in  Table  1.  The  parameters  were  derived  from  measurements 
taken  on-site.  Since  at  this  time,  full  three-dimensional  models  were  outside  the  capability  of  available 
computers,  two  approaches  to  the  problem  were  taken.  The  first  approach  was  to  compute  the 
waveguide  as  axially  symmetric.  The  second  approach  used  an  plane  wave  integral  transform  method 
which  assumes  invariance  in  one  spatial  dimension  of  the  waveguide.  In  this  case,  the  dimension  is 
perpendicular  to  the  domain  shown  in  Fig.  1.  These  solutions  were  compared  with  a  range 
independent  solution. 
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Figure  1:  The  domain  used  for  the  finite  element  model  of  TREX  Propagation  and 

Reverberation. 


Layer 

Water 

Mud 

Sand 

Compressional  Speed 

1525 

1514 

1769 

[m/s] 

Compressional  Att. 

Fran.-Garr. 

0.08 

0.3 

[dB/  ] 

Roughness  Spectrum 

Pierson-Moskowitz 

Flat 

Power-law 

(von  Karman) 

RMS  Roughness  [m] 

0.045 

0.0 

0.003 

3 

Density  [kg/m  ] 

1024 

1228 

2100 

Table  1.  Geoacoustic  parameters  used  in  the  finite  element  model  based  on  measurements  taken 

on  site. 


Analysis  of  Normal  Incidence  Bottom  Loss  Measurements  for  Range  Dependent  Geoacoustic 
Parameters 

Data  taken  at  the  Target  and  Reverberation  Experiment  in  May  2013  were  analyzed  and  compared  with 
bathymetric  and  high  frequency  backscatter  measurements  taken  on  site  by  deMoustier  and  Kraft. 
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Although  multiple  data  sets  were  taken,  the  bottom  loss  data  shown  in  Fig.  2  were  analyzed  for  initial 
bottom  loss  and  scattering.  Fig.  3  shows  the  match-filtered  data  as  a  function  of  depth,  which  was 
determined  from  the  sampling  frequency  and  measured  water  sound  speed  and  the  ping  number.  Each 
ping  had  an  associated  latitude  and  longitude  since  the  ROV  was  tracked  using  an  ultrashort  -baseline 
(USB)  acoustic  tracking  system.  This  allowed  the  ROV  data  to  be  directly  compared  to  the  multibeam 
survey  (MBS).  The  tracking  system  also  recorded  the  depth  of  the  ROV  and  with  the  bathymetry 
measured  from  the  MBS,  the  range  of  the  bottom  return  could  be  predicted.  This  is  shown  with  the 
heavy  black  line  in  Fig.  3.  The  bottom  loss  was  measured  by  considering  the  initial  bottom  return 
relative  to  a  calibrated  return  taken  in  the  laboratory  from  a  air/water  interface.  The  range  of  this  return 
was  determined  by  considering  the  maximum  of  the  signal  shown  with  the  blue  line.  Also,  considered 
was  the  scattering.  Its  relative  value  was  estimated  by  measuring  the  energy  in  the  signal  for  a  fixed 
time  after  the  initial  return.  The  fixed  time  interval  is  denoted  with  the  white  lines  in  Fig.  3.  These 
measurements  are  referred  to  as  “coda”  energy. 


Figure  2:  The  location  of  the  bottom  loss  measurements  relative  to  the  backscatter  of  the 

multibeam  survey. 
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Figure  3:  The  magnitude  of  the  match-filtered  data  of  each  ping  in  the  survey  shown  in  Fig.  2. 
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Measurements  of  Bottom  Loss,  Sediment  Structure  and  Interface  Roughness  at  the  Glider  Sensors  and 
payloads  for  Tactical  characterization  of  the  Environment  (GLISTEN  15)  Experiment 

ARL:UT  also  participated  in  the  Fliger  Sensors  and  payloads  for  Tactical  characterisation  of  the 
Environment  (GLISTEN  15)  experiment  conducted  in  August  and  September  of  2015  off  the  west 
coast  of  Italy.  Three  data  products  were  produced  by  ARL:  video  surveillance  of  the  ocean  bottom, 
normal-incidence  bottom  loss  data  which  is  also  an  indicator  of  scattering  and  layering  and  laser  line 
profiling  data.  The  expereimental  setup  is  shown  in  Fig.  4. 

The  data  collected  by  the  ARL:UT  system  will  aid  in  modeling  the  low-frequency  propagation  data 
taken  by  CMRE  as  part  of  the  experiment.  It  is  expected  that  ARL:UT  will  produce  a  finite  element 
model  of  the  propagation  to  understand  the  influence  of  range  and  depth  dependent  geoacoustic 
paramters.  To  those  aims,  data  was  taken  along  the  two  propagation  paths  of  the  experiment,  the 
north-south  track  which  was  relatively  range-independent  in  bathymetry  and  the  east-west  track  which 
was  range  dependent  in  bathymetry  with  more  than  a  200  m  depth  change.  The  locations  of  the  ROV 
deployments  are  shown  in  purple  relative  to  the  propagation  paths  in  Fig.  5.  In  the  figure,  north  is 
toward  the  top. 


Experimental  Setup 


ROV  Setup 


Deployment 


VESSEL  AT  ANCHOR  OR  IN  DP 


ROV  is  outfitted  with: 

-  Laser  profiling  system  for  interface 
roughness  measurement 

-  Normal  incidence  echo  sounder  from 
6-20  kHz 

-  ROV  camera 


Figure  4:  The  experimental  set-up  for  the  ARL:UT  ROV  for  the  GLISTEN  experiment. 
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Figure  5:  Locations  of  data  collection  for  the  ARL:UT  ROV  (shown  in  purple)  relative  to  the 

propagation  paths  for  the  GLISTEN  experiment. 


RESULTS 

Finite  Element  Modeling  for  Range  Dependent  Waveguides 

The  finite  element  propagation  model  was  compared  to  two  different  types  of  models  at  two  different 
frequencies  to  assess  the  effects  of  range  dependence  on  transmission  loss.  Shown  in  Fig.  6  is  the  axi- 
symmetric  FEM  compared  with  OASES  which  assumed  range  independent  parameters  and  simple 
cylindrical  spreading.  The  longitudinally  invariant  model  was  not  run  at  this  frequency  due  to  the 
large  number  of  out-of-plane  wavenumbers  required.  The  left  panel  shows  the  models  at  1  m 
increments  while  the  right  panel  shows  the  range-averaged  energy.  The  range-dependent  finite 
element  model  shows  significantly  more  loss  than  the  range-independent  case. 
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Figure  6:  Finite  element  propagation  loss  for  1900  Hz  over  domain  shown  in  Fig.  1.  In  the  right 
panel,  the  data  has  been  averaged  over  100  m  increments  to  show  the  change  in  level. 
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Propagation  loss  at  500  Hz  was  also  calculated  using  the  finite  element  model.  The  results  are  shown 
in  Fig.  7.  In  this  case,  the  longitudinally  invariant  model  was  also  calculated.  This  geometry  is  more 
suitable  to  the  physical  environment  since  the  sand  waves  were  relatively  longitudinal.  The  left  panel 
shows  the  data  averaged  over  200  m  increments.  In  this  case,  the  axi-symmetric  and  range- 
independent  cases  (as  shown  by  the  OASES  model)  are  largely  in  agreement.  However,  the 
longitudinally  invariant  finite  element  model  indicates  additional  loss  likely  due  to  out  of  plane 
scattering  by  the  sand  waves. 


Figure  7:  Propagation  model  for  500  Hz  over  domain  shown  in  Fig.  1.  The  longitudinally 
invariant  model  and  axi-symmetric  model  are  calculated  using  finite  elements.  In  the  right 
panel,  the  data  has  been  averaged  over  100  m  increments  to  show  the  change  in  level. 


Lastly,  low-frequency  model  of  reverberation  was  calculated  for  100  Hz  using  the  axi-symmetric  finite 
element  model  and  the  longitudinally  invariant  model.  The  results  are  shown  in  Fig.  8.  Although  the 
axi-symmetric  model  is  much  less  computationally  intensive  than  the  longitudinally  invariant  model,  it 
is  contaminated  by  backwards  propagating  cylindrical  waves.  In  the  case  of  propagation,  these  are 
much  lower  than  the  forward  going  energy  and  do  not  distort  the  result  as  drastically.  However,  for 
reverberation,  they  are  dominant  and  therefore,  the  axi-symmetric  model  is  unsuitable  for  reverberation 
modeling  in  which  the  geometry  is  longitudinal  such  as  the  case  with  sand  ripples. 
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Figure  8:  Finite  element  model  of  reverberation  calculated  at  a  center  frequency  of  500  Hz  over 

the  domain  shown  in  Fig.  1. 


Analysis  of  Normal  Incidence  Bottom  Loss  Measurements  for  Range  Dependent  Geoacoustic 
Parameters 

The  bottom  loss  and  coda  energy  data  derived  from  measurements  at  the  locations  shown  in  Fig.  2 
were  compared  with  the  depth  and  backscatter  from  the  MBS.  Shown  in  Fig.  9  is  the  bottom  loss  data 
compared  to  the  depth  (right  panel)  and  the  backscatter  (left  panel).  There  is  not  a  clear  correlation 
between  the  bottom  loss  which  at  these  frequencies  can  be  related  to  sediment  type  and  the  bathymetry 
or  backscatter.  This  may  have  implications  for  sediment  transport.  Conversely,  the  coda  energy  is 
compared  to  the  depth  (right  panel)  and  backscatter  (left  panel)  in  Fig.  10.  Here  a  clear  correlation  is 
seen  which  indicates  that  scattering  is  the  dominate  different  in  the  peaks  and  troughs  at  the  TREX  site 
and  that  scattering  is  also  the  main  mechanism  behind  the  large  variations  in  backscattering  from  the 
MBS  survey. 


Figure  9:  Bottom  loss  for  the  TREX  track  shown  in  Fig.  2  compare  with  MBS  values  of 
bathymetry  (right  panel)  and  backscatter  (left  panel). 
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Figure  10:  Coda  energy  for  the  TREX  track  shown  in  Fig.  2  compare  with  MBS  values  of 
bathymetry  (upper  panel)  and  backscatter  (lower  panel). 


Measurements  of  Bottom  Loss,  Sediment  Structure  and  Interface  Roughness  at  the  Glider  Sensors  and 
payloads  for  Tactical  characterization  of  the  Environment  (GLISTEN  15)  Experiment 

Since  the  GLISTEN  sea  test  was  so  recent,  only  preliminary  data  analysis  has  been  performed.  Shown 
in  Fig.  1 1  is  a  sample  of  the  normal  incidence  bottom  loss  data.  In  the  figure,  a  depth  dependent 
structure  is  evident.  The  sediment  was  a  mud  and  the  volume  inclusions,  likely  due  to  biology  based 
on  video  data,  are  clearly  visible.  Also,  a  lower  layer  is  visible  over  some  of  the  area.  These  data  have 
not  yet  been  accurately  assigned  latitude  and  longitude  based  on  a  USB  tracking  system. 
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Figure  11:  Preliminary  bottom  loss  data  from  the  GLISTEN  experiment. 


Shown  in  Fig.  12  is  an  example  of  the  laser  line  profiling  data.  Note  that  the  bottom  consisted  of  many 
burrows.  Brittle  stars  and  anemones  were  present  in  the  area  based  on  video  data.  These  laser  lines 
will  be  analyzed  to  provide  an  estimate  of  the  roughness  spectrum. 
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Data  from  Site  Near  S70 


Figure  12:  Laser  line  profiling  data  taken  near  site  S70. 
IMPACT/APPLICATIONS 

The  finite  element  reflection  loss  models  could  transition  into  a  new  high  frequency  and  low 
frequency  reflection  loss  (LFBL/HFBL)  data  curves  for  NAVO  based  on  site-specific  characteristics. 
It  is  also  a  benchmark  for  faster  approximate  methods  of  calculating  propagation  and  scattering.  The 
3D  LI  model  can  be  used  to  understand  propagation  and  reverberation  in  complex  environments.  An 
understanding  of  normal  incident  reflection  loss  is  critical  to  sediment  characterization  and  mine  burial 
prediction.  The  TREX13  measurements  will  serve  as  ground  truth  bottom  loss  and  interface  roughness 
measurements  for  reverberation  modeling.  GLISTEN  data  provide  insight  into  the  role  of  biology  on 
acoustic  propagation  and  scattering. 

RELATED  PROJECTS 

Under  the  iPUMA  and  SSAM  Sediment  Environmental  Estimation  (iSSEE)  program,  this  group  is  also 
developing  sediment  characterization  algorithms  for  AUV  sonars  based  on  the  measurements  and 
models  previously  developed  by  this  program.  Additionally,  the  models  developed  in  this  research 
will  be  used  to  increase  the  fidelity  of  sonar  trainers  under  the  High  Fidelity  Active  Sonar  Trainer 
(HiFAST)  program. 
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